The heavy-quark pole masses in the Hamiltonian 

approach 

^ . A.M. Badalian and A.I. Veselov, 

O I Institute of Theoretical and Experimental Physics, 

B.Cheremushkinskaya 25, 117218 Moscow, Russia 

00 , B.L.G. Bakker, 

Q"^ I Department of Physics, Vrije Universiteit, Amsterdam, The Netherlands 



(N 
> 
O 



Abstract 

From the fact that the nonperturbative self-energy contribution Cse 
to the heavy meson mass is small: Cse(&&) = 0; Cse(cc) = —40 MeV 
PP, strong restrictions on the pole masses mi, and rric are obtained. The 
analysis of the bb and the cc spectra with the use of relativistic (string) 
■ Hamiltonian gives m6(2-loop)= 4.78 ± 0.05 GeV and mc(2-loop)= 1.39 ± 

I 0.06 GeV which correspond to the MS running mass mb{rhb) = 4.19±0.04 

GeV and fhc{fhc) = 1.10±0.05 GeV. The masses uic and uob, which define 
^ l' the heavy quarkonia spin structure, are shown to be by ~ 200 MeV larger 

I , than the pole ones. 
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1 Introduction 

The spectrum of heavy quarkonia (HQ) is very rich and provides a unique oppor- 
tunity to study the static interaction in the infrared (IR) region, and hyperfine 
\ and fine structure effects. To use that opportunity one needs to know, besides 

such fundamental parameters as the string tension and the strong couphng, also 
the heavy-quark mass, which cannot directly be measured since a quark is not 
observed as a physical particle. Therefore the quark mass mg has to be deter- 
mined indirectly, e.g. from the study of hadronic properties like e^e~ bb, 
hadronic decays, and the QQ spectra. 

In the QCD Lagrangian the mass parameter depends on the renormalization 
scheme and by convention this current mass is taken in the MS scheme. In 
perturbation theory it is convenient to introduce the pole quark mass, i.e. the 
pole of the quark propagator, and at present the MS pole mass is known to 
three-loops OE]: 

mQ=mQ(mQ)<;i + ^^^+6|^) }, (1) 
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where the Lagrangian current masses, 



fhbifhb) = (4.25 ± 0.25) GeV, fh^{fhc) = (1.20 ± 0.20) GeV, (2) 



are known now with an accuracy of 17% (6%) respectively for the c quark {b 
quark). Most calculations of the pole masses nib and nic have been done in 
the QCD sum rules approach 4 , lattice QCD fSI, and different perturbative 
approaches |2IS1- 

For three decades many properties of HQ like the spectra, electromagnetic 
transitions, hadronic and semileptonic decays, were successfully studied in dif- 
ferent potential models (PM) 6 - 13 , however, the heavy quark masses used in 
PM are considered "to make sense in the limited context of a particular quark 
model" 12] i.e. as fitting parameters. 

However, in the last decade the situation has changed and in Ref. illj the 
relativistic (string) Hamiltonian was derived directly from the QCD Lagrangian, 
starting with the gauge-invariant meson Green's function in Fock-Feynman- 
Schwinger (FFS) representation. In Refs. mj-^SI it was established that for 
the orbital angular momentum L < 5 and not too large string corrections, as in 
HQ, the string Hamiltonian reduces to the well-known Hamiltonian Hq used in 
the relativized potential model (RPM) for many years |7| |H] : 



It follows from the derivation of Hq in Ref. that the mass in Q coincides 
with the MS running mass fhq{mq)^ if the perturbative interaction is neglected, 
or with the pole mass mq if the perturbative self-energy corrections are 
taken into account. Therefore this Hamiltonian can be used to extract the pole 
mass mq from the analysis of the HQ spectrum. 

Nevertheless, if one looks at the heavy-quark masses used in PM, a large 
variety of mt, and nic values can be found in different analyses: rric in the range 
1.30 GeV ^ 1.84 GeV and rrib in the range 4.20 5.17 GeV |Z|-|ini- The 
main reason behind this wide spread in the values of rn^ and rric (even for the 
same Hamiltonian _ffo) is the presence of a negative arbitrary constant Cq in 
the mass formula (or in the chosen static potential). We give three examples: 
in Ref. J nic = 1.327 GeV and Cq = is used by the Wisconsin Group; in 
Ref. ID nic = 1.628 GeV and Cq — —253 MeV (in both cases the Hamiltonian 
Eq. © was used); in Ref. ^ mc = 1.84 GeV and Co = -800 MeV are taken, 
i.e. the magnitude of Co is always larger for larger heavy-quark mass. 

The meaning of the constant Co was understood recently and in Ref. ^ 
it was shown that the negative contribution to the meson mass comes from 
the nonperturbative (NP) interaction of the color-magnetic moment of a quark 
(antiquark) with the background (vacuum) field. Moreover, this self-energy NP 
term Cse was analytically calculated with 10% accuracy |Tj (see Eq. (|59(l ^: 




(3) 



CsE(r?.i) 



4(7 



77(mg) 



(4) 
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for a quark and an antiquark with equal masses. In the expression Eq. Q 
is the pole mass which determines the factor rj{mq) Eq. (|57() while 

ujq{nL) = {^P^+^t},^^ (5) 

is the dynamical quark mass. For low- lying states in charmonium and bottomo- 
nium ujQ turns out to be ~ 200 MeV larger mg 

The essential fact (for light and heavy- light mesons) is that CsE,{nL) depends 
on the quantum numbers and just due to this the correct intercept of the Regge- 
trajectory has been obtained in Ref. ^Sl- In HQ the situation appears to be 
much more simple. The factor i]{'mq) in Eq. I^J depends on the flavor through 
the pole mass rriq and from the analytical expression (|57|l one obtains a small 
value: rjc = 0.35 0.27 for rric in the range 1.37 1.70 GeV and rji, = 0.07 for 
rUb = 4.7-;- 5.0 GeV. As a result Cse(&&) = —3 MeV, i.e. is compatible with 
zero, and Cse(cc) is also small: 

CsEibb) = 0; Cse(cc) = (-40 ± 10) MeV (6) 

Thus the self-energy contributions to HQ states are well defined and there- 
fore there is no opportunity anymore to vary the pole mass by introducing a 
fitting constant. We shall show in this paper that the condition ^ puts strong 
restrictions on the values of mbirric) needed to describe the bb(cc) spectrum. 
The extracted pole masses and rric in our analysis will be determined with 
an accuracy better than 60 MeV and the main uncertainty in their values comes 
not from the method used (for fixed string tension and the strong coupling (or 
Aqcd) the uncertainty is ±10 MeV) but from the uncertainty in our knowledge 
of the strong coupling in the IR region. We shall show that HQ spectra, in par- 
ticular high excitations and the recently discovered ID state in bottomonium, 
can give very important information about the strong coupling in the IR region. 

Our analysis of HQ spectra shows that in bottomonium mb(2-loop) < 4.70 
GeV and values mi, > 4.84 GeV turn out to be incompatible with the condition 
C'sE = 0. In charmonium the admittable rric values, TOc = 1.39 ± 0.06 GeV, 
appear to be rather small and agree with the one obtained by Narison with the 
use of the QCD sum rules for the (pseudo)-scalar current W . Our calculations 
of the HQ spectra are done with the use of only three fundamental quantities: 
the string tension, the QCD constant A(n/) and the pole mass toq. The main 
emphasis in our fit lies on the excited (not ground) states. 

This paper is organized as follows. In Section [21 the mass formula, following 
from the relativistic Hamiltonian, as well as the approximations to that, are 
presented and the notion of dynamical mass is introduced. In Section the 
static potential and the strong coupling in the IR region which is defined in 
background perturbation theory (BPT), are discussed. In Section 0] from the 
analysis of the bb spectrum (with special accent on high excitations) the pole 
mass mt, (2-loop) is obtained. In SectionElthe pole mass rric is extracted from the 
cc spectrum. In Section El our Conclusions are presented and in the Appendix 
the method and NP self-energy term are discussed. 
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2 The mass formula 



The string corrections are small in HQ and therefore the simplified form of the 
Hamiltonian Hji IT may be used (see Appendix): 

Hr = '—+LO+^+Vst.tic- (7) 
LO U) 

To derive this Hamiltonian in the FFS representation one needs to go over from 
the proper time r in the meson Green's function l|Sn|l to the actual time t and 
at this point a new variable uj{t) must be introduced 

This variable is the canonical coordinate and since Hji does not depend on its 
derivative, the requirement that the canonical momentum tt^^ — Q is preserved 
in time, corresponds to the extremum condition 

TTuj = {tt, Hr} = -q^^ 0. (9) 
From this extremum condition the operator lo, 



UJ^^p'^ + ml, (10) 

is defined as the kinetic energy operator. Substituting the definition (|10|l into 
the Hji |7J one arrives at the Hamiltonian Hq Eq. (PJ: 



Ho = 2.^p2+m2 + 1/3t,ti,(r), (11) 

which does explicitly not depend on the variable uj. However, to calculate 
different corrections to the meson mass, like spin and string corrections, or 
the self-energy corrections which can be considered as a perturbation and by 
derivation depend on the we shall use the approximation when for a given 
state the operator uj will be substituted by its average: 



This mass Lu{nL) can be called the dynamical mass since its difference with 
respect to the current (pole) mass mq is fully determined by the dynamics. 
Note that for vanishing pole mass the value of uj{nL) is finite and determines 
the constituent mass of a light quark. 

It is also important that perturbative corrections to the current mass, which 
are essential at small quark-antiquark separations, r < 0.1 fm |15[ I16| . are 
included in the pole mass mq occurring in Hq. On the other hand the static 
potential 14tatic('') is well defined at QQ separations r >Tg = 0.2 fm where Tg 
is the gluonic correlation length 1 . The eigenvalues of H^, denoted as Mo(nL), 

|2y^p2_^m| + y,tatic(r-)} V«L(r) = Mo{nL)^nL{r) (13) 
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together with the self-energy term Eq. JHJl define the heavy-meson masses. As 
shown in the Appendix, in bottomonium Cse — and therefore the spin- 
averaged mass M{nL) for a given bb state coincides with the eigenvalue Mo{nL): 

M{nL,bb) ^ Ma{nL), (14) 

while in charmonium from Eq. H57|l Cse — —40 MeV and 



M{nL, cc) = Ma{nL) + Cse- (15) 

There exist two approximations to the solution of the spinless Salpeter equation 
(SSE) leading to two approximations to the meson mass M{nL). First, the 
nonrelativistic (NR) approximation when the mass AI^-ix{nL) is given by: 

MNR(nL) = 2toq + Ef^lifUQ) + Cse, (16) 

where E^^{mQ) is the eigenvalue of the Schrodinger equation with the reduced 
mass equal to ^mg. 

There is also another, so called "einbein" approximation (EA) to the solu- 
tions of SSE H13|) . where the meson mass given by 

MEAinL) = u;{nL) + ^+ E^^lM + Cse (17) 

appears to be closer to the exact solution M{nL) than in the NR approxima- 
tion • In the binding energy i?„i(wQ) is defined as the solution of the 
Schrodinger equation with the reduced mass equal to ^uiQ{nL) (not the pole 
mass ^toq) while u!Q{nL) is to be defined from the selfconsistent equation: 

OMea „ rnl dE,,L{^) 

— = 0, or w„L = — 2- - w„L^ ■ (18) 

duo UnL OUJnL 

Through ujnL in EA the relativistic corrections are taken into account in the 
mass formula ()17|) . Moreover, owing to the special form of the mass formula 
H17|l the choice of uj{bb) = 5.0 GeV turns out to be compatible with the condition 
Cse = 0, while in NR approximation the admittable values of are about 200 
MeV smaller. 

It is worthwhile to notice that in bottomonium where both LijQ{nL) and mq 
are large, around 5 GeV, the diff'erence between NR, EA, and relativistic cases 
is small, \5r\ = M{nL) - M^^{nL) is about 10 ~ 20 MeV. In charmonium the 
difference depends on the quantum numbers and for high excitations can reach 
^ 100 MeV (see the discussion in Section 



3 Static potential 

The static potential contains perturbative and NP contributions where the NP 
linear potential can directly be derived from the meson Green's function if the 
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qq separation is larger than the gluonic correlation length. From the analysis 
of the Regge trajectories of light and heavy-light mesons the value of the string 
tension, a = 0.185 ±0.005 GeV^ is fixed while the perturbative interaction 
in coordinate space can be presentread in the form, 

Vp{r)^--^^^^^^, (19) 
3 r 

where astatic ('') is well known only in the perturbative region, i.e. at very small 
distances, r < 0.1 fm fT^IT?)] . 

However, the r.m.s. radii in bottomonium and charmonium, span a very 
wide range: 

i?(T(lS')) = 0.2 fm, i?(xb(lP)) 0.4 fm, i?(xfc(2P)) = 0.65 fm, 
R{T(4S)) = 0.9 fm, R{T{6S)) = 1.3 fm, (20) 

and 

RiJ/ip) = 0.4fm, R.{xc{lP)) = 0.6fm, R(^p{lD) = 0.8fm 
i?(V'(35)) = l.lfm, R{^j{4S) ^ 1.40fm. (21) 

Apparently, with the exception of T(IS') the sizes of these states lie outside the 
perturbative region. 

Therefore the problem arises how to define the strong coupling astatic (''') 
at all distances, in particular in the IR region. In PM it has always been 
assumed and later this fact has been supported by direct measurement of the 
static potential in lattice QCD, that the strong coupling freezes and reaches 
a critical (saturated) value at large r. Unfortunately, at present there is no 
consensus about the true value of acrit and different values were used. In the 
phenomenological analysis of Ref. |S]) acrit = 0.60 was used, but in analytical 
perturbation theory ^7] the large value acrit = 47r//3o = 1.4 appeared. In the 
background perturbation theory (BPT) which will be used here, acrit is smaller 
and fully defined by Aqcd TH^]. For the definition of astatic (j^) it is better 
to start with the vector coupling in momentum space: 

VB{q) = -AttCf^^ (22) 
9 

This background coupling a^iq) is defined in Euclidean momentum space at all 
q^, including q^ = 0, i.e. it has no Landau singularity. 



The logarithm 

tB=lr^^—rP^ (24) 



g2 + Mi 



A 
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contains the background mass Mb which appears due to the interaction of a 
gluon with the background field at small . This mass Mb = 1 GeV has the 
meaning of the lowest hybrid excitation: A/b — M{QQg) ~ M{QQ) and 
from the comparison with the static potential Mb determined on the lattice, 
was found to be equal to 1 GeV ^Hl- The ^b (|^4|I coincides in form with the 
parametrization of as{q) suggested in Refs.pH] where instead of the background 
mass Mb two gluonic masses {2mg) enter. However, the physical gluon cannot 
have a mass while Mb — M {QQgg) — M{QQg) is a well defined physical quantity 
(Mq is supposed to be large) and can be calculated in different theoretical 
approaches [2] and on the lattice. 

By definition aB('z) has the correct asymptotic freedom (AF) behavior at 
large and in this region the connection between the vector coupling asC?) 
and as((7) in the MS renormalization scheme is very simple, so that the QCD 
constant Ay (in the vector-scheme) can be expressed through Aj^ |22j: 

A^)=Aj-)expf^) (25) 



2/3, 



Here ai = ^ - ^n/. At present the value A^^l (2-loop) = 215 ± 15 MeV is 
established from high energy processes j^, while in quenched QCD the value 
Ai^ = 240 ± 20 MeV was calculated on the lattice El- Then from the relation 

MS ' * 

H25|) it follows that in quenched approximation the QCD constant in the vector 
scheme (|25|) has the value 

A^j"^ = 385 ± 30 MeV (26) 

and for nf ~ 5 

A^,^^ = 290 ± 30 MeV. (27) 

Our choice of Ab in this paper will be in accord with the numbers (|26|l and (|27|l . 

In coordinate space the background coupling (r) is defined as the Fourier 
transform of as {q) ■ 

oo 

, . 2 /■ sinqr 

aB[r) = - / dq aB[q), (28) 

71- 7 q 



so that the perturbative part of the static potential is 

V^ir) = -t^ (29) 
and the static potential is the sum of Vb (r) and the NP linear potential 

V^.tatic(r)=ar-^^. (30) 
3 r 

In phenomenology the approximation where the coupling aB{r) is constant is 
often used. This approximation is valid, because at distances r > 0.4 fm the 
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Table 1: The effective coupling acs{nL) (|33|l for different bb and cc states for 
the static potential (jSOI) with nib = 4.78 GeV, = 1.45 GeV, ctcrit = 0M7;a = 
0.185 GeV2, A(4) = 360 GeV (n/ = 4). 



State 


IS 2S 3S 4S 5S 6S 


aos{nL, bb) 
Ucff {nL, cc) 


0.386 0.419 0.427 0.430 0.431 0.432 
0.441 0.447 0.446 0.445 


State 


IP 2P 3P ID 2D IF 


acsinL, bb) 
acsinL, cc) 


0.474 0.480 0.482 0.508 0.508 0.520 
0.510 0.505 - 0.531 0.526 



coupling as(r) is already approaching the freezing, or critical, value. Also at 
large distances, r > 1.2 fm, the confining linear potential, due to pair creation, 
is becoming more flat |24ll25j and this effect can be important not only for light 
mesons but also for high excitations in charmonium. 

Note that the critical values of aB{q) and asir) coincide, 

"Big = 0) = Q!B(r oo) = acrit, (31) 

and their characteristic values in two- loop approximation are given below, (A'^"^ ^ = 
A^)) 

385 MeV) = 0.428 
370 MeV) = 0.510 

340 Mev) = 0.515. (32) 

Our calculations show that the bottomonium spectrum appears to be rather 
sensitive to the AF behavior while in charmonium the approximation ae = 
const can be used with good accuracy. It is also instructive to look at the 
effective coupling aofi for different cc and bb states, which can be defined as 

{^)^^ = -MnL){r-X^ (33) 

and is dependent on the quantum numbers (see Table and about 20% is 
smaller than acrit- 

From Table^one can see that in bottomonium due to the different character 
of the wave functions the effective coupling is smaller for the nS states and rather 
large for orbital excitations like the nD states. 



4 Bottomonium 

To extract the pole mass m?, the bb spectrum will be studied here as a whole. 
We mostly ignore the ground state - the T(IS') mass, for which high perturba- 
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Table 2: The spin-averaged masses M{nL) (in GeV) in bottomonium for the 
potential H30|l with the parameters H35|) (NR case). 



State 


Set I 


Set II 


Set 3 


experiment 


IS 


9.460 


9.406") 


9.379"^) 


9.4603± 0.0003 


2S 


10.013 


10.001a) 


9.988 


10.233± 0.0003 


3S 


10367 


10.359 


10.359 


10.3552±0.0005 


IP 


9.900 


9.900 


9.900 


9.901± 0.0006 


2P 


10.267 


10.267 


10.270 


10.260± 0.0006 


ID 


10.150 


10.156 


10.161 


10.1612± 0.0016 




Above 


BB threshold, A/th= 


10.558 GeV 


4S 


10.659 


10.647 


10.649 


10.5800± 0.0035 


5S 


10.917 


10.900 


10.902 


10.865±0.008 


6S 


11.146 


11.131 


11.132 


10.019±0.008 



°) The mass M{1S) increases by an amount of ~ 50 80 MeV if the AF 
correction is taken into account. 

tive corrections can be important [3], but rather concentrate on the following 
experimental splittings pill4| 

M^ogilD) - M{1P) = Af (l^Ds) - M{1P) = 261.8 ± 1.8 MeV (exp) 
M{2P) - Af (IP) = 360.0 ± 1.2 MeV (exp). (34) 

It is important that all IP, 2P, and ID states lie below the BB threshold and 
have no hadronic shifts. Also for these splittings the small relativistic corrections 
are partly cancelled and the calculations can be done either with the use of SSE 
or in NR approximation. 

First, we consider the case with astatic — const and give the hh spectrum in 
Table 121 for three sets of parameters with different mfc: 

Set I mb= 4.727 GeV, ct = 0.20 GeV^ astatic = 0.3345, 

Set II mb= 4.765 GeV, ct = 0.19 GeV^ astatic = 0.390, 

Set III mb= 4.778 GeV, ct = 0.188 GeV^ astatic = 0.415. (35) 

In all cases Cse = is taken in the mass formula H14|) . 
From the masses presented in Table |5] one can see that 

1. For small m?, — 4.727 GeV (Set I) the mass Af(l_D) appears to be about 
10 MeV lower than the experimental value even for very large a = 0.20 
GeV^. 

2. For Set II and Set III almost identical fits are obtained with exception of 
the ID state where good agreement with experiment can be reached only 
for a larger value of the coupling, as for Set III. 
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Table 3: The dynamical mass 0Jb{nL) (in GeV) and the difference between the 
dynamical mass and the current quark mass (in MeV) for SSE with Cornell 
potential (nib = 4.78 GeV, a ^ 0.185 GeV^, astatic = 0.4125). 



State 


IS 


2S 


3S 


4S 


5S 


6S 


IP 


2P 


ID 


uj(nL) 


5.043 


5.008 


5.028 


5.057 


5.088 


5.127 


4.959 


4.989 


4.962 


ujinL) — nib 


263 


228 


248 


277 


308 


347 


179 


209 


182 



3. The spin-averaged ID — IP splitting, 

A = A/cog(li?) - Mcog(lF), (36) 

has remarkable properties - it is practically independent of the relativistic 
correction Sr and small variations of the string tension, and therefore A 
can be considered the best and very stable criterium to determine the 
critical value of the strong coupling as well as the pole mass nib- 

4. In NR approximation Eq. (|16|l and for SSE for the 6-quark masses mb < 
4.70 GeV or nib > 4.8 GeV the condition Cse = cannot be combined 
with a reasonably good fit to the bb spectrum. 

However, if one uses EA Eq. H17() instead of the NR mass formula Eq. H16|l 
the values of the dynamical mass ujb{nL) are larger and the difference 
u!b{nL) — irib varies in the range (180 300 MeV) (see Table 3), from 
which one can see that the dynamical mass uJb{nL) is slightly different for 
different nL states. The bb spectrum calculated with uJb{nL) = 5.0 GeV 
and rub = 4.78 GeV with the use of the mass formula l(T7j) gives values 
coinciding within ±5 MeV with those from Table 2. 

Thus from our fits, when the coupling is taken constant, the extracted value 
of the pole mass 

rub = 4.76 ± 0.02 GeV (astatic = const) (37) 

is obtained. 

The picture does not change much if the AF behavior of aB{r) in two-loop 
approximation is taken into account. 

However, in this case the admissible values of nib appear to be larger by ~50 
MeV. The bb spectrum for nib = 4.82 GeV and nib — 4.83 GeV for the number 
of the flavors nf = 4, 5 and also in quenched approximation is presented in Table 
4. 

As seen from Table 0] good agreement for Sets A and B is obtained for the 
IP and ID states but for the IS* level the mass is ~ 10 MeV higher than the 
M(T(15)) value. This fact can be connected with the contribution of the 3-loop 
perturbative correction which is neglected here. 
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Table 4: The bb spectrum defined by the mass formula (|16|l for the spinless 
Salpeter equation. 





Set A 


Set B 


Set C 






rub = 4.816 


nib = 4.83 GeV 


rub = 4.817 






a = 0.185 GcV2 


cr = 0.19 GeV2 


cr 0.185 GeV^ 




State 


A(5) = 360 MeV 


A(4) =390 MeV 


A(o)(l-loop)=365 MeV 


experiment") 


IS 


9.471 


9.478 


9.470 


9.460 


2S 


10.025 


10.032 


10.023 


10.023 


3S 


10.376 


10.386 


10.375 


10.355 


IP 


9.900 


9.901 


9.900 


9.900 


2P 


10.271 


10.278 


10.266 


10.260 


ID 


10.159 


10.162 


10.152 


10.161 



) The experimental errors of the masses are given in Tab. [3 



Our conclusion is that for the vector constant A^^^ = 390 MeV, or A*^^^ < 365 
MeV, and a = 0.186 ± 0.004 GeV*, the extracted pole mass lie in the narrow 
range 

m6(2 - loop) = 4.81 ±0.02 GeV. (38) 

Thus, if the AF behavior of q:b (r) is taken into account the extracted pole mass 
is about 50 MeV larger than m;, for astatic = const. 

Then combining Eqs. H37|l and H36|l for different choices of ae one obtains 
that the extracted pole mass of the b quark lies in the range 

mb = 4.78 ± 0.05 GeV. (39) 

Then by definition of the two- loop pole mass {uf — 5) where in the relation 
Q the parameter ^2 is 

^2(n/ - 5) = |-1-0414^ (^1 - + 13.4434| , (40) 

and the sum over k extends over the flavors Qk which are lighter than Q, 
one finds ^2{nf = 5) = 9.6-;- 9.7, and for as{fhb) = 0.217 it follows from |^ that 

rfibifhb) = (4.19 ± 0.04) GeV, (41) 

This number for the MS mass H41|l appears to be in good agreement with the 
conventional value from Ref. but has smaller theoretical error, ^ 50 MeV, 
than the number quoted in Eq. 

5 Charmonium 

The cc spectrum has several differences in comparison to bottomonium. 
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Table 5: The cc spin- averaged masses M(nL) (in MeV) in R and NR cases 
with the same static potential (astatic = 0.42; a = 0.18 GeV^), rric = 1-41 GeV 
and different Cse (MeV). The quantity 6r is the relativisitic shift M{nL) — 
M^^inL). 



State 


Experiment 


R case 


NR case 


Sr 






Cgg = — 35 


CgT = -57 




IS 


3067± 0.7 


3067 


3067 (fit) 





2S 


3673± 6 


3661 


3688 


-27 


IP 


3525± 0.6 


3528 


3510 


18 




Above DD threshold 


ID 


3770± 2.5 


3823 


3812 


12 




3872±1.2 








2P 




3965 


3984 


-19 


IF 




4067 


4067 





3S 


4040±10 


4082 


4141 


-59 


2D 


4159±20 


4200 


4227 


-27 


4S 


4415±6 


4433 


4527 


-94 



First, the self-energy contribution to the mass M{nL) is nonzero, about —40 
MeV, being practically the same for different nL states, and therefore can be 
taken constant for all states with the accuracy 1-^3 MeV. 

Secondly, relativistic corrections arc not small in charmonium and M{nL) 
is in the relativisitc case (R) always smaller, so 

SninL) = M{nL) - Af^^(nL) (42) 

is negative. Note that the self-energy term must be the same in both cases . 
However, if the spin-averaged mass of the IS state, Af(lS')=3067 MeV, is used 
for the fit, then Cse are different in R and NR cases and 5R{nL) has irregular 
behavior (see TablelSj). From TableEJ where Cse — —35 MeV in the R case and 
Cs^^ = —57 MeV in the NR case, one can see that 5R{nL) is positive (~ 10-;- 20 
MeV) for the IP and ID states; equal to zero for the IF state, and negative for 
the nS states, 2P and 2D, and higher states. It is important that for the 45' 
(35*) state |(5(ni)| is large, ~ 100 MeV (60 MeV) and therefore the cc spectrum 
has to be calculated with a relativistic Hamiltonian. 

The interesting observation is that while the ground state mass is fitted well, 
the relativistic corrections to M-^iiinL) for the IP and ID states turn out to 
be positive, (since a negative value for in NR with larger magnitude was 
taken) . 

A third difference refers to the choice of a-B{r). Since the cc states have 
larger sizes than the hh ones the AF behavior of ae {r) appears to be less impor- 
tant in charmonium and the approximation Q:B(^)=constant is valid with good 
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Table 6: The spin-averaged masses in charmonium in R case for the static 
potential with the parameters given in Eq. (|43|l . 



State 


Set A 


Set B 


Experiment 


IS 


3067 


3067 


3067± 0.7 


2S 


3660 


3668 


3673± 8 


IP 


3528 


3510 


3525± 0.6 


Above DD threshold 


ID 


3823 


3805 


3871.8± 1.2 








3770.± 2.5 


2D 


4199 


4198 


4159 ± 20 


3S 


4080 


4109 


4040± 10 


4S 


4424 


4459 


4415± 6 


2P 


3964 


3954 





accuracy. For example for two sets of parameters: 

Set A mc = 1.42GcV, ct = 0.18 GcV, q^b = 0.42; CgE = -35 MeV 

Set B mc = 1.42GeV, ct = 0.185 GeV. A^^) = 360 MeV, Cse -30 MeV 

(43) 

close values of M(nL) in the relativistic case are obtained (see Table EJ. 

As seen from TableElthe higher excitations, like the 3S, 4S, and 2D states, 
lie ~ 40 MeV higher than the experimental values. All these states have large 
r.m.s. radii: R{3S) = 1.1 fm, R{AS) ^ 1.4 fm, and R{2D) = 1.4 fm. At 
such distances the confining potential is flattening due to quark-antiquark pair 
creation and it results in a correlated shift of the radial excitations down as 
it takes place for light mesons This phenomenon can be illustrated taking 
instead of the linear potential agr the modified confining potential a{r) r which 
was proposed in Ref. [25) . 

t \ ti \\ fi \ exp(^(r-a)) 

o-(r) ^CTo(l-7o/('')) with /(r) = — — ^ (44) 

B + exp(Y^ {r - a)) 

with the parameters 

0- = 0.185 GeV^ 70 = 0.40, a = 6.0 GcV"\ B = 20. (45) 

For this set of parameters the cc spectrum (R case) is given in Table [7| together 
with the one for the standard linear potential o-Qr with ctq = 0.185 GeV^. The 
value astatic = 0.42 is taken in both cases. 

From Table [71 one can see that for the modified potential a(r)r the masses 
M(4S') and MiZS) of the radial excitations are shifted down by 50 MeV and 
20 MeV respectively, and turn out to be closer to the experimental values. It is 
also worthwhile to look at the dynamical masses for the low-lying states which 
are larger than the pole mass by the constant amount 

ujc{nL) -mc = 220 ~ 250 MeV, (46) 
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Table 7: The comparison of the spin-averaged masses in charmonium (R case) 
for confining aor potential and modified potential H44() (toc — 1-42 GeV, Cse = 
-42 MeV, astatic = 0.42 in both cases). 



State 


(To = const = 


a — a{r) with 


experiment 




= 0.185 GeV2 


parameters Eq. 1451) 




IS 


3068 


3067 


3067 


2S 


3670 


3664 


3672 


3S 


4097 


4077 


4040±10 


4S 


4454 


4403 


4415±6 


IP 


3535 


3530 


3525±0.6 


2P 


3979 


3965 




ID 


3835 


3828 


3779±25 


2D 


4217 


4194 


4159± 20 



Table 8: The dynamical masses ujc{nL) for the potential with the parameters 
Eq. (gSJ and = 1.42 GeV (Set A). 



State 


IS 


2S 


3S 


4S 


IP 


2P 


ID 


2D 


WcijiL) in GeV 


1.65 


1.69 


1.74 


1.76 


1.63 


1.69 


1.66 


1.77 



while for high excitations this difference can reach 300 340 MeV (see Table 
inj. The observed difference between the dynamical and the pole mass can be 
essential for such physical characteristics as the hyperfine and fine-structure 
splittings, which are determined by the dynamical mass [21] and it causes a 
decrease of the hyperfine splitting, e.g. for the 25* state in charmonium |27| . 

In our analysis the best fit to the cc spectrum together with the correct 
choice of the self-energy contribution Eq. (O gives the pole mass rric in the 
range 

mc(2-loop) = 1.39 ± O.OlGeV (theory) ± 0.04(aB). (47) 

Then the MS running mass Eq. jSJ («/ = 4, the coefficient ^2 — 10.5) from 
Eq. (gZI) is 

TOc(mc) = (1.10 ± 0.05)GeV. (48) 

The value obtained turns out to be in good agreement with the conventional 
value for rhc{fri(.) (2j, but has a smaller theoretical error. 

6 Conclusion 

Our study of the bb and cc spectra has been performed with the use of the rel- 
ativistic Hamiltonian Hq and correct NP self-energy contribution to the meson 
mass. 
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By derivation the kinetic part of Hq contains the pole quark mass mg and 
it can directly be extracted from the analysis of the QQ spectrum. In our study 
all meson masses are expressed through only two parameters: the string tension 
and the QCD constant A^"^^ (in the Vector scheme). In charmonium the strong 
coupling aBir) can be approximated by a constant with good accuracy. 

The spin-averaged splittings like ID-IP and 2P-1P in bottomonium and 2S- 
IP and IP-IS in charmonium appear to be very sensitive to the chosen freezing 
(critical) value of the strong couphng. A good description of the HQ spec- 
tra was reached only if ctcrit was taken rather large, acrit — 0.55 ± 0.02 while 
the constant value for acs taken in the Coulomb potential is about 20% smaller. 

From our analysis one can conclude that 

(i) The dynamical quark mass tjf,(wc) is about 200 MeV larger than the pole 
mass mb(mc) for low-lying states. This difference should be taken into account 
when the spin structure in heavy quarkonia is studied, and it is especially im- 
portant in charmonium. 

(ii) The pole masses, TO6(2-loop)=4.78 ± 0.05 GeV and = 1-39 ± 0.06 
GeV, were extracted from our fit to the QQ spectra which correspond to the 
MS running masses: fhb{fhb) = 4.19±0.04 GeV and ffidTfic) = 1.10±0.05 GeV. 
The numbers obtained are in good agreement with the conventional values but 
have smaller theoretical error. The error we found is small, because in our 
analysis only one parameter, A (or acrit), is actually varied while a second 
parameter-the string tension-was taken the same as for light mesons. 

A Relativistic Hamiltonian 

Here we present the principal steps to derive the Hamiltonian Hr Eq. (jjj) and 
the NP self-energy term Eq. taken from Refs. ^E]- The starting point 
is the gauge-invariant meson Green's function written in FFS representation 
1211 with the use of the QCD action 



In Eq. H49|) the averaging goes over the background field and Gq{x, y){Gq{x, y)) 
is the Euclidean quark (antiquark) Green's function 



GM{x,y) = {Tr Pi Gq{x,y)T2 Gq{x,y))B 



(49) 



Gq{x,y) 



{rhq + D) 



-1 



{niq - DUml - D')-}y 



(50) 



OO 








where the factors T^a, "T^b, and IZp are given by, 



■a 
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Pf exp ^ J gaf_t^Ff^„dT^ 



(51) 



(The factors corresponding to the antiquarks are defined similarly.) Here Pa, Pb, 
and Pp are the ordering operators of the matrices a^, i?^, and F^^ respectively, 
and 



SE aH 



o'liu = (l/4*)(7/i7i' ~ Ivlij) represents the interaction of the quark (antiquark) 
color magnetic moment with the NP field strength F^^. 

In Eq. H50() the kinetic energy term K is defined as the integral over the 
proper time r: 



K 



rrig s - 



+ 7 (if^)dr. 



(53) 



The quark moving along the trajectory z^{t) interacts with the field of the 
valence gluon and by its color charge with the NP background field _B^. 

In Eq. H53|l the quantity rriq is the Lagrangian current mass usually taken 
in the MS renormalization scheme. The factor TZa is responsible for the stan- 
dard perturbative corrections to the quark mass rhq (as in Eq.(^), i.e. for the 
appearance of the pole mass in the QCD Action (Hamiltonian) 3 . Finally, the 
factors TZa and TZa (from the quark and the antiquark) provide the perturbative 
static interaction |25) . 

The other two factors TZb and TZb (from the quark and the antiquark) in 
GM{x,y) Eq. ifSUIl are responsible for the full NP (string) dynamics and were 
considered in detail in Ref. ^I], where after several steps the meson Green's 
function was presented in following form 



GM{f) = J du! dv dr exp{—Aji), 



(54) 



where the action Ar in coordinate space is expressed through two auxiliary fields 
Lu and ly. Since this action (see Ref. ^J) does not depend on the derivatives 
uj and v, the integration over uj, v in Eq. H54|l is equivalent to the canonical 
quantization of the Hamiltonian Hr which corresponds to the action Ar. It 
results in the following Hamiltonian, 



Hr 



UJ 



L2 

,2 



df3 Piy{f3) 



- 2 2 

-cr r 



d(3 



dpyifi), 



(55) 
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where the field operator uj{t) is defined as uj{t) — and t is the actual time. 
With the use of the extremal conditions {vq = <^f) and considering the string 
corrections as a perturbation (the procedure is described in Ref. [T3j) one obtains 
the Hamiltonian Hq. The terms TZp and TZp (from quark and antiquark) provide 
the NP self-energy contribution Cse (gauge-invariant) to the meson mass 0, 
where for the quark (antiquark) the self-energy correction Ato^ to the pole mass 
TUq appears to be expressed through the string tension a and the factor rj: 

ATOg(m„) ^ -—r]{m„) (56) 

TT 

The factor r]{mq) was calculated in analytical form in Ref. and for rriq > Tg, 
where Tg is the gluonic correlation length {5 = T~^), r]{mq) is given by the 
expression 



3my 5\2ml + 5^) 

'^(™.) = - (^2 _ ^2)5/2 --^t-- 5 + („,2 _ ^2)2 (57) 

A straightforward calculation gives ri{mB = 5.0) = 0.07; rj{nnq — 1.70) — 0.24 
and 'q{mq = 1.40) = 0.30. Then the unperturbed part of the string Hamiltonian 
Eq. {Tj) acquires the correction Eq. H56II 

7^2 ^2 _j_ /\y^2 

i^o ^ T +ojq + Fstatic ^Ho + Cse (58) 

LOq 

with the self-energy correction 

Ami 4:a 
Cse = — ^ = %{mq). 59) 

UJq TTUJq 

If this self-energy correction is considered as a perturbation, the operator cD in 
Eq. ifS^ can be replaced by the average of this operator Eq. i.e. by the 
dynamical mass 

'p2+m2)„i (60) 

and through cuq the NP self-energy term Cse appears to be dependent on the 
quantum numbers nL. However, in bottomonium 

C'sEibb) ^ -3 MeV (61) 

is small and can be neglected in the mass formula. In charmonium however, 
for rric = 1.40 GeV the factor ?]c = 0.29 and the value of Cse = -40 MeV is 
obtained which is practically the same for different nL states because of the 
weak dependence of LudnL) on the quantum numbers (see Table 
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